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Atmospheric neutrino results from Super-Kamiokande and Kamiokande 
— Evidence for oscillations — 
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New atmospheric neutrino results from Super-Kamiokande are presented. Results from Kamiokande on upward 
going muons are also presented. All these data, together with the Kamiokande atmospheric neutrino data give 
evidence for neutrino oscillations. Two flavor <-> v T oscillations, with large sin 2 2# and Am 2 in the region of 
10~ 3 to 10~ 2 eV 2 , explain all these data. 
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1. Introduction 

Cosmic ray interactions in the atmosphere pro- 
duce neutrinos. The prediction of the abso- 
lute flux has an uncertainty of ±20%. However, 
the flavor ratio of the atmospheric neutrino flux, 
(V M + Vfj)/{v e + v e ), has been calculated to an ac- 
curacy of better than 5% in a broad energy range 
from 0.1 GeV to higher than 10 GeV. The cal- 
culated flux ratio has a value of about 2 for en- 
ergies <lGeV and increases with increasing neu- 
trino energy. For neutrino energies higher than 
a few GeV, the fluxes of upward and downward 
going neutrinos are expected to be nearly equal; 
the geomagnetic field effects on atmospheric neu- 
trinos in this energy range are expected to be 
small because the primary cosmic rays that pro- 
duce these neutrinos have rigidities exceeding the 
geomagnetic cutoff rigidity (~ 10 GeV/Ze). 

The Vfi/ve ratio has been measured in deep un- 
derground experiments by observing final-state 
lcptons produced by charged-current interac- 
tions of neutrinos on nuclei, v + N — > I + 
X. The measurements are reported as R = 
(p/e)data/(p/e)MC, where (i and e are the num- 
ber of muon-like (/i-like) and electron-like (e- 
like) events observed in the detector for both 
data and Monte Carlo (MC) simulation. This 
ratio largely cancels experimental and theoret- 
ical uncertainties, especially the uncertainty in 
the absolute flux. The observed values of R 
by Kamiokande § were significantly smaller 



than unity for both below and above 1 GeV en- 
ergy ranges. Consistent results were obtained by 
IMB-3(<1.5GeV) @ and Soudan-2 @. 

Neutrino oscillations have been suggested to 
explain the small values of R. For a two-neutrino 
oscillation hypothesis, the probability for a neu- 
trino produced in a flavor state a to be observed 
in a flavor state b after traveling a distance L 
through a vacuum is: 

'1.27Am 2 (eV 2 )£(km 



P a ^ b = sin 2 29 sin 2 



(1) 



Eu(GeV) 

where E v is the neutrino energy, 6 is the mixing 
angle between the flavor eigenstates and the mass 
eigenstates, and Am 2 is the mass-squared differ- 
ence of the neutrino mass eigenstates. For detec- 
tors near the surface of the Earth, the neutrino 
flight distance, and thus the oscillation probabil- 
ity, is a function of the zenith angle of the neu- 
trino direction. Vertically downward-going neu- 
trinos travel about 15 km while vertically upward- 
going neutrinos travel about 13,000 km before in- 
teracting in the detector. The broad energy spec- 
trum from a few hundred MeV to about 100 GeV 
and this range of neutrino flight distances makes 
measurements of atmospheric neutrinos sensitive 
to neutrino oscillations with Am 2 down to 10 -4 
eV 2 . 

The zenith angle dependence of R measured 
by the Kamiokande experiment]^] at high ener- 
gies, together with the small R values, has been 
cited as evidence for neutrino oscillations. Based 



2 



on these measurements, Kamiokande obtained 
the allowed parameter regions of neutrino oscil- 
lations. Because of the relatively small statistics, 
both <-> v e and <-> v T oscillations were al- 
lowed. 

Recently, a long baseline reactor experiment, 
CHOOZ excluded the <-> v e solution of the 
atmospheric neutrino problem. 

The upward-going muons observed in under- 
ground detectors are the products of neutrino in- 
teractions in the rock. The mean neutrino energy 
of these events are of the order of 100(10) GeV for 
through-going(stopping) events. These events are 
used for the independent check of the neutrino os- 
cillation analysis of the lower energy atmospheric 
neutrino data. 

We present the analyses of atmospheric neu- 
trino events from Super-Kamiokande. In addi- 
tion to the events whose vertex positions are in 
the fiducial volume of the detector, we present 
the upward-going muon results from Kamiokande 
and Super-Kamiokande. We observed small val- 
ues of R and a zenith angle dependent deficit of 
/i-like events. While no combination of known 
uncertainties in the experimental measurement 
or prediction of atmospheric neutrino fluxes is 
able to explain the data, a two-neutrino oscilla- 
tion model of <-> v x , where v x may be v T or 
a new, non-interacting "sterile" neutrino, is con- 
sistent with the observed R values and zenith an- 
gle distributions. These data and the neutrino 
oscillation interpretation were further supported 
by a small (upward stopping muons) / (upward 
through-going muons) ratio and zenith-angle de- 
pendent deficit of upward through-going muon 
flux. From these measurements, we conclude that 
the atmospheric neutrino data give evidence for 
neutrino oscillations. 

2. Super-Kamiokande detector 

Super-Kamiokande is a cylindrical 50 kton wa- 
ter Cherenkov detector located at a depth of 2700 
meters water equivalent in the Kamioka Observa- 
tory in Japan. The detector consists of an inner 
detector surrounded by an outer detector on all 
sides. The inner and outer detectors are opti- 
cally separated by a pair of opaque sheets. 11146 



50 cm <fi photomultiplier tubes (PMTs), instru- 
mented in all surfaces of the inner detector, de- 
tect Cherenkov photons radiated by relativistic 
charged particles. 1885 20 cm <fi PMTs are instru- 
mented in the outer detector. The outer detector 
is useful for identifying entering cosmic-ray muons 
and measuring exiting particles produced by neu- 
trino interactions occurring in the inner detector. 
Pulse height and timing information from each 
PMT are recorded and used in the data analysis. 
The trigger threshold for electrons is 5.7 MeV/c 
at 50% efficiency. 

For a description of the Kamiokande detector, 
see Ref.||. 

3. Fully and partially contained events 

Super-Kamiokande observed a total of 4353 
fully-contained (FC) events and 301 partially- 
contained (PC) events in a 33.0 kiloton-year ex- 
posure. FC events deposit all of their Cherenkov 
light in the inner detector while PC events have 
exiting tracks which deposit some Cherenkov 
light in the outer detector. For the present anal- 
yses, the neutrino interaction vertex was required 
to have been reconstructed within the 22.5 kilo- 
ton fiducial volume, defined to be > 2 m from 
the PMT wall. The number of FC+PC events 
observed so far in Super-Kamiokande was about 
4 times larger than that in Kamiokande. 

FC events were separated into those with a sin- 
gle visible Cherenkov ring and those with multiple 
Cherenkov rings. For the analysis of FC events, 
only single-ring events were used. Single-ring 
events were identified as e-like or /i-like based on 
a likelihood analysis of light detected around the 
Cherenkov cone. The FC events were separated 
into "sub-GeV" {E ms < 1330 MeV) and "multi- 
GeV" (E ms > 1330 MeV) samples, where E ms is 
defined to be the energy of an electron that would 
produce the observed amount of Cherenkov light. 
E ms = 1330 MeV corresponds to -1400 MeV/c 
for muons. 

In a full-detector Monte Carlo simulation, 88% 
(96%) of the sub-GeV e-like (/i-like) events were 
v e (y^) charged-current (CC) interactions and 
84% (99%) of the multi-GeV e-like (/i-like) events 
were v e (v^) CC interactions. PC events were es- 
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Data 


Monte Carlo 


sub-GeV 






single-ring 


2389 


2622.6 


e-like 


1231 


1049.1 


/x-like 


1158 


1573.6 


multi-ring 


911 


980.7 


total 


3300 


3603.3 


R = 0.63 ± Qm(stat.) ± 0.05(sys.) 


multi-GeV 






FC events 






single-ring 


520 


531.7 


e-like 


290 


236.0 


[i-like 


230 


295.7 


multi-ring 


533 


560.1 


total 


1053 


1091.8 



PC events 



total(=^-like) 301 371.6 

R FC+PC = 0.65 ± Q.05(stat.) ± 0.08(sys.) 



Table 1 

Summary of the sub-GeV, multi-GeV and PC 
event samples observed in 33 kiloton-year expo- 
sure of the Super-Kamiokande detector. The data 
are compared with the Monte Carlo prediction 
based on the neutrino flux calculation of Ref . p] . 



target, 3% from the neutral current cross sec- 
tions, 2% from particle identification, 1% from 
the absolute energy calibration, 0.6% from the 
vertex fit and fiducial volume cut, less than 
0.5% from non-neutrino background events and 
1.5% from statistical uncertainty in the Monte 
Carlo. Adding these errors in quadrature, the 
total systematic uncertainty was 8%. The sys- 
tematic uncertainty of R for the multi-GeV sam- 
ple, obtained in a similar way, was 12%. Ta- 
ble | summarizes the measured R values in 
Super-Kamiokande. These results are consis- 
tent with the Kamiokande R values, which were 
0-60 tom( stat -) ± 0.05(sys.) for the sub-GeV 
data and 0.57 ±?$(stat.) ± 0. 07 (sys.) for the 
multi-GeV datag. 

The Super-Kamiokande data have been ana- 
lyzed independently by two groups, making the 
possibility of significant biases in data selection 
or event reconstruction algorithms remote[Q,||]. 
Given the statistical error in R especially for the 
sub-GeV sample, statistical fluctuation can no 
longer explain the deviation of R from unity. As- 
suming that the systematic error has the gaus- 
sian form, we estimate the probability that the 
observed fi/e ratios could be due to statistical 
fluctuation is less than 0.001% and less than 1% 
for sub- and multi-GeV samples, respectively. 



timated to be 98% CC interactions; hence, all 
PC events were classified as /^-like, and no single- 
ring and particle identification requirements were 
made. Table [l] summarizes the number of ob- 
served events for both data and Monte Carlo as 
well as the R values for the sub-GeV and multi- 
GeV samples. Further details of the detector, 
data selection and event reconstruction used in 
this analysis are given in Ref. |t],|8| . 

Flavor ratio 

Super — Kamiokande measured significantly 
small values of R in both the sub-GeV and multi- 
GeV samples. Several sources of systematic un- 
certainties in these measurements were consid- 
ered(0,|). For the sub-GeV sample, they were: 
5% from uncertainty in the predicted v^/Ve flux 
ratio, 3.5% from the CC neutrino interaction 
cross sections and nuclear effects in the H2O 



Zenith angle distribution 

The /x-like data from Super-Kamiokande ex- 
hibit a strong up-down asymmetry in zenith angle 
(O) while no significant asymmetry is observed in 
the e-like data, see Figure [l]. For further analyses, 
we define up-down ratio U / D where U is the num- 
ber of upward-going events (—1 < cos 9 < —0.2) 
and D is the number of downward-going events 
(0.2 < cos 8 < 1). The ratio is expected to be 
near unity independent of flux model for E v > 1 
GeV, above which effects due to the Earth's mag- 
netic field on cosmic rays are small. Based on a 
comparison of results from our Monte Carlo sim- 
ulation using different flux models as inputs, 
treatment of geomagnetic effects results in an un- 
certainty of roughly ±0.02 ~ 0.03 in the expected 
U/D values for e-like and ^t-like sub- GeV events 
and less than ±0.02 for multi-GeV events. These 
two flux calculations do not assume the existence 
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Super-Kamiokandc 



Kamiokande 



MC 



Data 



Data 



e-like 

Sub-GcV, <400MeV/c 
Sub-GcV, >400McV/c 
Multi-GeV 
/it-like 

Sub-GeV, <400MeV/c 
Sub-GcV, >400MeV/c 
Multi-GeV (FC+PC) 



1.00±0.04±0.03 
1.02±0.04±0.03 
1.01±0.06±0.03 

1.05±0.03±0.02 
1.00±0.03±0.02 
0.98±0.03±0.02 



1.20±g;io =•= °- 03 
1.10tg;iJ±0.03 
0.93±g;i| ± 0.02 



1.03±g;iJ ± 0.02 
0.65±g;$j ±0.01 
0.54±g : g| ± 0.01 



1.29] 



1+0-27 
-0.22 

76+ - 22 

u - ' °-0.18 

1 OO+0.39 
1 -°°-0.30 



i 1O+0.31 
1.18_ 24 

i-09±g:S 



Table 2 

Summary of the up/down ratio, [//I?, for e-like and /z-like events from Super-Kamiokande and 
Kamiokande. The systematic errors of the Kamiokande experiment are not shown, but are similar to 
those of Super-Kamiokande. 
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Figure 1. Zenith angle distributions observed in Super-Kamiokande for; (a)sub-GeV e-like, (b)sub-GeV 
/z-like, (c)multi-GeV e-like and (d) multi-GeV (FC+PC) /i-like events. CosG —1 means down-going 
particles. The histograms with shaded error bars show the MC prediction with their statistical errors for 
the no neutrino oscillation case. The dotted histograms shows the Monte Carlo prediction for «-> v T 
oscillations with sin 2 20 =1 and Am 2 =2.2xl0~ 3 eV 2 . 
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of 1 km mountain over the Super-Kamiokande 
detector. The rock reduces the neutrino flux be- 
cause muons are stopped before they can decay in 
flight. We estimated the effect of the presence of 
rock on the flux. U/D changed about 2% for the 
multi-GeV events. The change was much smaller 
for the sub-GeV events. 

Studies of decay electrons from stopping muons 
show at most a ±0.6% gain difference, i.e., mea- 
sured energy difference, between up-going and 
down-going particles. This gain difference caused 
0.9%, and 0.7% uncertainty in U/D for the multi- 
GeV e-like and /i-like events. This uncertainty is 
much smaller for the sub-GeV sample. A con- 
tamination of non-neutrino background such as 
down-going cosmic ray muons could have direc- 
tional correlation. The maximum contribution 
to the uncertainty in U/D from the contamina- 
tion was estimated to be ±1.0%, ±0.1%, ±0.5% 
and ±2.0% for the sub-GeV e-like, /i-like, multi- 
GeV e-like and /i-like events, respectively. From 
these studies, the total systematic uncertainties 
in U/D for the data and MC are summarized 
in Table [|. In the present data, U/D for e-like 
events is consistent with expectations. U/D for 
high momentum /i-like events significantly devi- 
ates from unity, while U / D for low momentum 
/i-like events is consistent with unity. The av- 
erage angle between the final state lepton direc- 
tion and the incoming neutrino direction is 55° 
at p = 400 MeV/c and 20° at 1.5 GeV/c. At 
the momentum range below 400 MeV/c, the pos- 
sible up-down asymmetry of the neutrino flux is 
largely washed out. We have found no detector 
bias differentiating e-like and /i-like events that 
could explain an asymmetry in /i-like events but 
not in e-like events ||. 

The U/D value for the multi-GeV FC+PC 
/i-like events, 0.54 tplos ^ 0-01 deviates from 
unity by more than 6 standard deviations. This 
value is also consistent with the Kamiokandc U / D 
value for the multi-GeV FC±PC /i-like events, 



0.58 



+0.13 
-0.12- 



These numbers, which are close to 



0.5, suggest a near maximal neutrino mixing. 

Neutrino oscillation analysis 

We have examined the hypotheses of two-flavor 
Vu. <-> v e and Va <-» v T oscillation models using 



a x 2 comparison of the Super-Kamiokande data 
and Monte Carlo, allowing all important Monte 
Carlo parameters to vary weighted by their ex- 
pected uncertainties |ll|. The data were binned 
by particle type, momentum, and cos 8. A x 2 is 
defined as: 



x 2 = E 

cosO,p 



N Da ta - N MC {sin 2 29, Am 2 , e 3 ) 



E 



(2) 



where the sum is over five bins equally spaced in 
cos and seven momentum bins for both e-like 
events and /i-like plus PC events (70 bins total). 
The statistical error, a, accounts for both data 
statistics and the weighted Monte Carlo statis- 
tics. Noata is the measured number of events in 
each bin. Nmc (sin 2 29, Am 2 , ej) is the expected 
number of Monte Carlo events and is a function 
of sin 2 29, Am 2 and tj. tj are parameters which 
are related to the systematic uncertainties. The 
parameters (and their uncertainties) considered 
in this analysis are: overall normalization (25%, 
but this was fitted as a free parameter), E v spec- 
tral index (0.05), sub-GeV R (8%), multi-GeV 
R (12%), relative normalization of PC to FC 
(8%), L/E v (15%), sub-GeV (2.4%) and multi- 
GeV (2.7%) up-down ratios. See Ref.|I| for more 
details. 

For <-> v e , effects of matter on neutrino 
propagation through the Earth were included 
following Ref. JO]. Due to the small number 
of events expected from r-production (15 to 20 
events were expected in the present FC±PC sam- 
ple), the effects of r appearance and decay were 



neglected in simulations of 



A global 



scan was made on a (sin 2 29, log Am 2 ) grid mini- 
mizing x 2 with respect to uncertainty parameters, 
ej, at each point. 

The best-fit to <-» v T oscillations, 

X 2 mm{phys) = 65.2/67 DOF, was obtained at 
(sin 2 29, Am 2 ) = (1.0, 2.2 x 10~ 3 eV 2 ) 
inside the physical region (0 < sin 2 29 < 
1). The best-fit values of the Monte Carlo 
uncertainty parameters were all within their 
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expected errors for this point. The global 
minimum occurred slightly outside the phys- 
ical region at (sin 2 29, Am 2 ) = (1.05,2.2 x 
10- 3 eV%xLn(un P H V s) = 64.8/67 DOF. The 
90% C.L. allowed region is located at X 2 min ( phys ) + 
5.0, based on the minimum inside the physical 
region [Q. The allowed region is shown in Fig- 
ure H In the region near x 2 minimum, the \ 2 
distribution is rather flat and has many local min- 
ima so that inside the 90% interval the best-fit 
Am 2 is not well constrained. The \ 2 increases 
rapidly outside of the 90% C.L. region. We ob- 
tained x 2 — 135/69 DOF, when calculated at 
sin 2 29 — 0, Am 2 = (i.e. assuming no oscil- 
lations). 

The Kamiokande allowed regionQ obtained by 
the contained event analysis is also shown in 
Figure [| The Super-Kamiokande allowed re- 
gion favors lower Am 2 than that of Kamiokande. 
However, the allowed regions from both experi- 
ments have a region of overlap. One of the rea- 
sons for the difference in the allowed region is 
due to the zenith angle distribution of the sub- 
GeV /z-like events. The U/D values of the /i-like 
events from Kamiokande and Super-Kamiokande 
are summarized in Table ^. It should be noted 
that Kamiokande U / D value was consistent with 
unity for /z-like events in the momentum range 
of sub-GeV, >400MeV/c but Super-Kamiokande 
observed smaller U/D value for the same sam- 
ple; a 2.5cr difference in U/D for this sub-sample. 
This difference results in a difference in the fa- 
vored Am 2 region in these two experiments, since 
the energy region observing the small U/D is di- 
rectly related to the determination of Am 2 , 

For the test of «-> v e oscillations, the 
Super-Kamiokande data resulted in a rela- 
tively poor fit; xL„ = 87.8/67 DOF, at 
(sin 2 2(9, Am 2 ) = (0.93, 3.2 x 10~ 3 eV 2 ). The ex- 
pected U/D value of the multi-GeV e-like events 
for the best-fit «-> v e oscillation hypothesis, 
1.52, differs from the measured value, 0.93 ig 12 ^ 
0.02, by 3.4 standard deviations. We conclude 
that the v p <-> v e hypothesis is not favored. 

The zenith angle distributions for the sub- 
and multi-GeV samples are shown in Figure |l|. 
The data are compared to the Monte Carlo ex- 




Mass^MeV/c 2 ) 

Figure 2. Invariant mass distribution observed 
in Super-Kamiokande for two ring events with 
"e-like, e-like" particle-identification results and 
with no /i-decay signal. The histogram with 
shaded error bars show the Monte Carlo predic- 
tion with their statistical errors. Cut region of 
7r°-mass is also shown. 



pectation (no oscillations, hatched region) and 
the best-fit expectation for v p «-> v r oscillations 
(dashed line). The oscillated Monte Carlo well 
reproduces the zenith angle distributions of the 
data. 

In Super-Kamiokande, the oscillation param- 
eters were also estimated by considering the R 
measurement and the zenith angle shape sepa- 
rately. The 90% C.L. allowed regions for each 
case overlapped at 1 x 10~ 3 < Am 2 < 4 x 10~ 3 
cV 2 for sin 2 26 = 1. 

As another check of the Super-Kamiokande 
data, we studied 7r° events. For «-> v T os- 
cillations, the number of neutral-current (NC) 
events should be unchanged by neutrino oscil- 
lations. As a NC sample, we selected "7r°"-like 
events. The selection criteria were: 2 ring events, 
both rings should be e-like, no /x-decay signal 
and the invariant mass calculated from the charge 
and direction of the two rings assuming two 7's 
should be between 90 and 190MeV/c 2 . We ob- 
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served 210 7r°-like events. See Figure g, which 
shows the invariant mass distribution for the two 
ring events. A clear excess at ~140MeV/c 2 is 
seen. The fraction of NC events in this sam- 
ple is estimated to be 82%. For «-> v r os- 
cillations, the number of CC v e events (<~ e-like 
events) should also be unchanged by neutrino os- 
cillations. Therefore, the (7r°/e) ratio of the data 
should agree with the same ratio of the Monte 
Carlo without oscillations for the «-> v r case. 
We obtained: (tt /e) data /(ir° /e) MC — 0.93 ± 
0.07(stat) ± 0.19(syst) (Preliminary). The re- 
sult is consistent with the <-> v T interpretation 
of the data. 

4. Upward going muons 

Energetic atmospheric f M 's passing the Earth 
interact with rock surrounding the detector and 
produce muons via CC interactions. Because 
of the atmospheric muon background, it is dif- 
ficult to select neutrino induced downward going 
muons. On the contrary, upward going muons are 
essentially neutrino origin. 

Upward going muons can be categorized into 
two types. One is "upward through-going muons" 
which are the events which enter into the detec- 
tor and exit, and the other is "upward stopping 
muons" which enter the detector and stop in the 
detector. 

Upward through-going muons 

The mean energy of neutrinos which produce 
upward through-going muons is about lOOGeV. 
Kamiokande observed 372 upward going muons 
during 2456 detector live days. The selection cri- 
teria were; cosO = — 1 ~ —0.04 and the minimum 
track length in the inner detector of 7 meters [fbif. 
The minimum (mean) energy loss of these muons 
in the inner detector is 1.6(3.0)GeV. The average 
detection efficiency was 97%. With the require- 
ment, cos© < —0.04, the background contamina- 
tion was negligible. The observed flux of upward 
going muons was 1.94±0.10(stai.)tQ;og(sys.) x 
10~ 13 cm~ 2 sec~ 1 sr~ 1 . The expected flux 
based on the calculated flux of Ref. was 
2.46±0.54(i/ieo.) x lO-^cm^sec^sr" 1 . Fig- 
ure shows the zenith angle distribution of 
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Figure 3. Zenith angle distribution of 

upward through-going muon flux observed in 
Kamiokande. Inner(outer) error bars show statis- 
tical (statistical+uncorrelated experimental sys- 
tematic) errors. The solid histogram shows 
the expected flux for the null neutrino oscilla- 
tion case. The dashed histogram shows the ex- 
pected flux for the <-> v T oscillation case with 
sin 2 26> =1.0, Am 2 =3.2xl0~ 3 eV 2 and a =1.00. 



the upward through-going muon flux observed in 
Kamiokande. 

In Super-Kamiokande, 617 upward through- 
going muon events were observed during 537 de- 
tector live daysjlj]. The selection criteria were; 
cosO <0, two outer detector cluster correspond- 
ing to the muon entrance and exit points and 
track length of a muon inside the inner detector 
should be longer than 7 m. The minimum (mean) 
energy loss of these muons in the inner detector is 
1.6 (6) GeV. The average detection efficiency of 
these events was estimated to be >99%. The va- 
lidity of this efficiency was tested using the real 
down-going cosmic-ray muons by assuming the 
up-down symmetry of the detector. 

The number of background events, 4.6, was es- 
timated by extrapolating the zenith-angle distri- 
bution of cosmic ray muons with cos© =0^0.08. 
The background events are expected only in the 
cos© = —0.1 ~ bin and are subtracted in 
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Figure 4. Zenith angle distribution of up- 
ward through-going muon flux observed in Super- 
Kamiokande. Error bars show statistical + un- 
corrected experimental systematic errors. Esti- 
mated background is subtracted. The solid his- 
togram shows the expected flux for the null neu- 
trino oscillation case based on the calculated flux 
of Ref||. The dashed histogram shows the ex- 
pected flux for the <-> v T oscillation case with 
sin 2 20 =1.0, Am 2 =2.5xl(T 3 eV 2 and a =1.12. 



x 10 13 cm 2 sec 1 sr x ). Fig- 



further analyses. The observed flux of upward 
going muons was 1.75±0.07(sfat.)±0.09(sys.) x 
10 _ cm~ 2 sec~ 1 sr~ 1 . The expected flux 
based on the calculated flux of Ref . |p[ ( |Io|] ) 
was 1.88±0.42(t/ieo.) x 
(2.01±0.45(i/ieo. 

ure |^ shows the zenith-angle distribution of the 
upward through-going muon flux observed in 
Super-Kamiokande. 

Because of the difference in the mean track 
length of muons, the observed and predicted 
fluxes are different between the two experiments. 
However, the observed/predicted flux ratios are 
consistent between the two experiments within 
the measurement errors. The measured zenith- 
angle distributions in these two experiments have 
similar shape; both experiment observed lower 
flux near the vertical direction compared with 



the predicted flux distributions. The x 2 values 
of the comparison of the shape of the zenith an- 
gle distributions of the data and predictions are 
21.3/9 DOF and 18.7/9 DOF for the Kamiokande 
and the Super-Kamiokande data, respectively. 

Analyses of neutrino oscillations were carried 
out in these experiments. To test the oscillation 
hypothesis, a \ 2 is defined as: 



X 



10 

E 

COS& 



(f>Data ~ a ■ 4>mc (sin 2 28, Am 2 ) 



a — 1 



(3) 



where a is the statistical and systematic error 
in the observed flux (ipData) and ot and a a are 
an absolute normalization factor of the expected 
flux(0Mc( sm2 28, Am 2 )) and its uncertainty, re- 
spectively. a a was taken to be ±22% by adding 
the theoretical uncertainty and correlated exper- 
imental errors. The uncertainty in the absolute 
neutrino flux (±20%) was the dominant source of 
<7 Q . A minimum x 2 wa s calculated by changing 
a for each (sin 2 2(9, Am 2 ). Since the contained 
data prefer <-> v T oscillations and since the 
CHOOS experiment § excluded the <-> v e os- 
cillation parameter region relevant to the atmo- 
spheric neutrino data, only <-> v r oscillations 
were tested. 

The Kamiokande data had 

xLn at (sin 2 2(9, Am 2 ) = (1.0,3.2 x 10~ 3 eV 2 ) 
with a =1.00. The xLin value was 12 - 8 / 8 D0F - 
Xmin f° r the Super-Kamiokande data occurred 
at (sin 2 2<9, Am 2 ) = (1.00,2.5 x 10" 3 eV 2 ) with 
a =1.12. The xL™ value was 7 - 3 / 8 D0F - Thc 
allowed regions did not change significantly for 
two assumptions of the neutrino flux]9]|,[Tc|], sug- 
gesting that the zenith angle and energy depen- 
dences of the neutrino flux at high energies are 
understood well. Figure ^ shows the allowed re- 
gions obtained from these upward going muon 
data. The allowed regions are larger than those 
of contained events. But they overlap the allowed 
regions obtained by the contained data. 

Upward stopping muons 

Because of the large detector dimension of 
Super-Kamiokande, a substantial fraction of up- 
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ward going muons stop in the detector. The mean 
neutrino energy of the upward stopping muons is 
about lOGeV, which is substantially lower than 
that of upward through-going muons. There- 
fore, in some neutrino oscillation parameters, it 
is expected that an observed (stopping/through- 
going) flux ratio of upward going muons is differ- 
ent from the calculated ratio. 

Super-Kamiokande observed 137 upward stop- 
ping muons during 537 detector live days Jl5| . 
The selection criteria were similar to those of 
upward through-going muons except for a re- 
quirement of one outer detector cluster corre- 
sponding to the entrance point. The detection 
efficiency was estimated to be 99%. The esti- 
mated number of cosmic-ray background events 
was 13.2. It was estimated that these back- 
ground events were mostly in the cosO = 
—0.1 ~ bin. The background events were 
subtracted from this bin for further analyses. 
We then calculated the (stopping/through-going) 
flux ratio (= 5ft). The observed value was 
5ft =0.22±0.023(siat.)±0.014(sys.), while the pre- 
dicted value was 5ft =0.39±0.05(t/ieo.). The ob- 
served value was substantially smaller than the 
prediction. Figure || shows the zenith angle dis- 
tribution of the (stopping/through-going) flux ra- 
tio. Also shown in the same figure are the pre- 
dicted ratios with and without neutrino oscilla- 
tions. Clearly, the predicted distribution with 
neutrino oscillations agrees well with the data 
within the systematic uncertainty (±14%). 

The 90% C.L. allowed region of the neutrino 
oscillations was estimated by a x 2 test: 



E 

COS& 



^Data -P -Kmc (sin 2 29, Am 2 



(3-1 

<T0 



(4) 



where a is the experimental (mostly statistical) 
error in the observed ratio (Sftuata) and and 
<rp are a uncertainty factor of the expected ratio 
(5ftMc(sin 2 29, Am 2 )) and its 1 a error, respec- 
tively. a a was taken to be ±14% by adding the 
theoretical uncertainty and correlated experimen- 
tal systematic errors. 




Figure 5. Zenith angle distribution of the 
(stopping/through-going) ratio of the upward go- 
ing muon flux observed in Super-Kamiokande. 
Error bars show statistical + uncorrelated exper- 
imental systematic errors. The solid histogram 
shows the expected ratio for the null oscilla- 
tion case. The dashed histogram shows the ex- 
pected ratio for the <-> v T oscillation case with 
sin 2 26> =1.0 and Am 2 = 3xl0~ 3 eV 2 . The ex- 
pected ratio has ±14% correlated uncertainty. 



Xmin occurred at (sin 2 20, Am 2 ) = (1.0,3.4 x 
10" 3 eV 2 ). The xL„ value was 1.3/3 DOF. Fig- 
ure |^ shows the allowed region obtained from this 
analysis. Again, the allowed region is larger than 
those of contained events. But it overlaps the al- 
lowed regions obtained by the contained and up- 
ward through-going muon data. 



5. Conclusions 

Both the zenith angle distribution of /i-like 
events and the (/x/e) values observed in Super- 
Kamiokande were significantly different from the 
predictions in the absence of neutrino oscilla- 
tions. While uncertainties in the flux predic- 
tion, cross sections, and experimental biases are 
ruled out as explanations of the observations, the 
present data are in good agreement with two- 



flavor 



or v u 



oscillations. This 
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Figure 6. The allowed neutrino oscillation 
parameter regions obtained by Kamiokande and 
Super-Kamiokande. Inside of each curve is al- 
lowed at 90% C.L.. The (l)thick-black and 
(2)thick-gray curves show the regions obtained by 
contained event analyses from Super-Kamiokande 
and Kamiokande, respectively. The ^black- 
dotted and (4)gray-dotted curves show the re- 
gions obtained by upward through-going muons 
from Super-Kamiokande and Kamiokande, re- 
spectively. The (5)thin-black curve shows the re- 
gion obtained by the (stopping/trough-going) ra- 
tio analysis of upward going muons from Super- 
Kamiokande. 



conclusion is consistent with the results from 
Kamiokande on contained data analysis and sup- 
ported by the upward-going muon results from 
Super-Kamiokande and Kamiokande. Two ex- 
periments, Super-Kamiokande and Kamiokande, 
give consistent data and various techniques point 
to a common parameter region of neutrino oscil- 
lations: Am 2 should be in the range around 10 
~ 10" 2 eV 2 and sin 2 26>>0.8. We conclude that 
the atmospheric neutrino data, especially from 
Super-Kamiokande, give evidence for neutrino os- 
cillations. 

The Super-Kamiokande experiment is sup- 
ported by the Japanese Ministry of Education, 



Science, Sports and Culture and the United Sta- 
tus Department of Energy. The Kamiokande ex- 
periment was supported by the Japanese Ministry 
of Education, Science, Sports and Culture. 
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